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ABSTRACT 

Using simple arguments we show that the canonical thin keplerian accretion 
disks cannot smoothly match any plain advection dominated flow (ADAF) model. 
By 'plain' ADAF model we mean the ones with zero cooling. The existence of 
sonic points in exact solutions is critical and imposes constraints that cannot be 
surpassed adopting 'reasonable' physical conditions at the hypothetical match 
point. Only the occurrence of new critical physical phenomena may produce a 
transition. We propose that exact advection models are a class of solutions which 
don't necessarily involve the standard thin cool disks and suggest a different 
scenario in which good ADAF solutions could eventually occur. 

Subject headings: : accretion disks — black hole physics — hydrodynamics 

1. Introduction 

Recently, the need to explain black hole candidate spectra and the low luminosity 
of some galactic nuclei, supposed to contain a black hole, prompted more attention 
to advective accretion disk models (Narayan and Yi 1995) The basic idea is that the 
gravitational energy content doesn't go into local radiation emission as in standard thin 
and cool keplerian disks (let us call SSD) (Shakura and Sunyaev 1973), but goes into 
heating and into radial motion. Therefore a new class of disks has been proposed: ADAF 
(Advection Dominated Accretion Disks). These disk models take into account the advection 
and pressure terms in the fluid-dynamics equations. Indeed, exact solutions of accretion 
flows (including viscosity, but only for isothermal cases) have been found years ago by S. 
Chakrabarti (1990). Literature has been enriched by many kinds of such advective disks 
with different physical properties, see Bisnovatyi-Kogan (1999) for an extensive review. 
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We will examine the very simple ADAF models initially proposed (Narayan and Yi 1995; 
Narayan, Kato and Honma, 1997; Igumeshev, Abramowicz and Novikov, 1998; Dullemond 
and Turolla,1998; Lu, Gu and Yuan, 1999), i.e. the ones without any cooling. As stated by 
the authors, although this model may appear very crude, it contains the basic ingredients: 
the role of pressure gradient and of the advection term. Analytical solutions have been 
obtained only for the self similar case (Narayan and Yi 1995), but this approximation is 
not tenable for the case of BH accretion, since all standard ADAF solutions rely on the 
conditions at the sonic point, which is close to the BH. 

Other solutions have been typically obtained with numerical integration of the 
steady state equations. The most used technique is the relaxation method with boundary 
conditions given by 'reasonable' physical approximations. In all works it is under laying the 
meaningful but 'a priori' idea that necessarily the ADAF must spring out of a canonical 
disk. So all the numerical procedures compel the solutions to fit the standard kepler disk. 
It is worthwhile to add that it is well known that the relaxations methods can relax onto 
unphysical solutions. 

A simple energy consideration gives a cutting argument against the possible smooth 
convergence of ADAF to SSD disks. The fact that self similar solution by Narayan and 
Yi cannot be automatically connected to cold thin disks has already been indicated by 
Kato and Nakamura (1998). Here we point out that even the more sophisticated (not self 
similar) ADAF cannot have an asymptotic approach to the SSD models. We conclude that, 
possibly, standard 'plain' ADAF models are 'numerically forced' solutions and the claimed 
asymptotic approach to a keplerian disk is not demonstrated. 
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2. Basic equations 



In the 'plain' ADAF model it is assumed that it doesn't contain any diffusive term 
and the cooling processes are irrelevant, as occurring on timescales much longer than the 
advection fall time. No heat conduction is considered as in Narayan and Yi (1995), Narayan, 
Kato and Honma (1997), Dullemond and Turolla (1998). Then the basic equations are the 
ones given by classical viscous fluid dynamics. We adopt newtonian physics and the black 
hole force on the fluid can be derived from the appropriate Paczyhski and Wiita potential 
(1980). This is also a standard procedure in many papers. On the other hand the classical 
approach makes more clear the basic physical point without inessential complications due 
to an exact relativistic treatment. 

Therefore an exact disk solution has to obey the following equations: 

Mass conservation 



M— rYiV r = const 



Radial momentum equation 



dv r 
r dr 



E dr 




where £ and II are the vertically integrated density and pressure, 



Stress definition 



T r4> = T l r ~Q r ~ 



where r/ = apaZdisk with a 




is the sound speed, 



Angular momentum equation 



pv r dX 
r dr 



1 d 

r 2 dr 
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that leads to 

M A = 4irr 2 T r(f) + J 
Vertical equilibrium gives the disk thickness 



7 - — 

^disk 



Energy equation 



de _ de PI d(rv r ) $ 
dt r dr p r dr p 



with e thermal energy per unit mass, $ the dissipation function given by : $ = 



2 

TV, 



Equation of state 

P=(l-l)pe 

Qk is the keplerian angular velocity, all other terms have their standard meaning 
accepted in the referenced papers. 



3. Boundary conditions 

With the help of the energy equation, that will be discussed in the next section, our 
problem is reduced to solve a set of 3 algebraic equations (vertical equilibrium, mass, energy) 
and two differential equations (radial momentum and angular momentum). The existence 
of sonic points, requiring the continuity of the variables and of the speed derivatives, links 
these boundary conditions. In particular, once the sonic point is determined, then also 
the derivative of the radial speed at that point is fixed. The derivation of the analytical 
formula of at the sonic point is algebraically complicated but straighforward. The usual 
assumption |^ = at the sonic point is reasonable only for sonic point close to the BH 
and for low viscosities, and it is not true in general. Due to the different approaches to the 
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numerical solutions (relaxation from fixed inner and outer points, standard Runge-Kutta 
integration from sonic point) there is no consensus on the exactness of the adopted boundary 
conditions. Our solutions are computed integrating the equations starting from the sonic 
point with the appropriate analytical derivative values calculated by the De L'Hopital rule. 
Artemova et al. (2001) using a similar treatment of the sonic point conditions have obtained 
subkeplerian flows starting from canonical disks. However their results don't disprove 
our conclusions, since they refer to disks obtained taking into account radiation transfer 
and using different viscosity prescriptions. Apart any further consideration on boundary 
conditions, we show that some crucial conditions are imposed by the energy equation. 



4. The energy argument 

An interesting point, apparently not exploited in ADAF context, although appearing 
in many papers ( Hoshi 1984, Honma 1996) is the fact that it is possible to have a constant 
energy property, reformulating the energy equation in the floowing way. 



div 



l>\ ( ^v 2 • M'ir) • // ) - v : T 







h — — + £ is the enthalpy, T is the stress tensor. 

This expression, after insertion of the stress definition, integrating over Z, and using 
the mass conservation equation, in general gives the following relationship, even in the 
differential form of the viscosity tensor: 

1 2 M* a 2 A 2 AA e const 

A e is the angular momentum at the inner edge of the disk, usually set equal to the 
value corresponding to the last stable orbit around the Black Hole. 
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The energy equation is valid over the whole ADAF solution up to the eventual crossing 
point with the keplerian disk. A smooth transition between ADAF and SSD disk requires a 
continuous smooth limit of the ADAF quantities to the SSD quantities: temperature, radial 
speed, density, etc.; 

What is crucial in this formula is that, at the crossing point of the solutions, where the 
angular momentum is keplerian, the following relation is obviously valid 

2r 2 r 2 

since it is reasonable and common to choose A e <C A fc . Now, if the ADAF solution has 
to smoothly match a standard keplerian disk at a radial distance far from the BH (r ^> r g ) , 
we can substitute A with A&; another reasonable and usual assumption is that at the match 
point the flow be subsonic (Igumenshchev et al. 1998). 

Therefore at the crossing point r + we must have : 

« 2 =K + 2Mi + 2L = K+ ^ 



(7 — 1) r + 2r + 2 2 r + 

Now, requiring that the energy constant be of the same (or of the same order) of the 
energy per unit mass of a standard keplerian disk, i.e. K = —\^f^- , we have : 

a 2 GM* 



(7 - 1) r + 

corresponds to the virial temperature , much larger than the standard temperature 
of a canonical thin keplerian disk. We have a strong inconsistency: the ADAF temperature, 
at the matching point, is always much larger than the canonical keplerian disk temperature 
(Frank, King and Raine 1992). 
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Indeed Lu et al. (1999) use our same numerical approach (i.e. starting from sonic 
point) and apparently they find cases with a smooth convergence from ADAF to SSD; but 
it is worth nothing that their parameter space, corresponding to Narayan ADAF model, has 
zero measure. Those configurations lie on a line in their parameter space R s , A e . Therefore 
they can be obtained only for a single A e if the sonic point is given. 

5. Solutions 

We performed several integrations using the direct method, starting from the sonic 
point inwards and outwards. Figure 1. shows the result of the Runge-Kutta integration 
of our equations. It is plotted the angular momentum per unit mass versus the radial 
distance for a =0.005, 0.01, 0.02. As the viscosity increases the angular momentum crosses 
the keplerian angular momentum at distances closer to the BH. Figure 2. shows the sound 
speed and the Mach number for the same a values. The sound speed at the crossing point 
doesn't go to extremely low values. The values we obtain are exactly the ones predicted by 
the energy equation and go to a finite limit at a finite radius located beyond the keplerian 
crossing point. The angular momentum of ADAF solutions crosses the kepler angular 
momentum value in a non smooth way. The exact solution continues beyond that point 
and stops at a larger distance with disk conditions (temperature, density, radial speed) not 
corresponding to the canonical keplerian disk. Our correct ADAF solutions are 'per se' 
standing solutions. This kind of solutions is a particular case of the more general solutions 
drawn by Chakrabarti and Titarchuck (1995) and Chakrabarti (1996), who explored a 
more general case including an 'ad hoc' cooling treatment. We add a further comment on 
the shock or no shock solutions, since many authors consider the Chakrabarti solutions 
untenable as the shock he finds had not been confirmed by their calculations. As has been 
pointed out in early works (Liang and Thompson, 1980) there are many sonic points in 
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the flow solutions. The flow on a black hole has no rigid surface and shock occurs only in 
solutions passing through the outer sonic point as clearly demonstrated by Chakrabarti 
(1990) and Chakrabarti and Molteni (1993). All the ADAF solutions we are discussing 
pass through the inner sonic point, so it is quite obvious that no shock can occur in these 
branches: the flow is supersonic only after the sonic point which is close to the black 
hole horizon. Lanzafame et al. (1998) tested these solutions with a time dependent code 
and found that these solutions are stable, while the standard ADAF solutions have been 
demonstrated unstable by numerical simulations (Igumenshchev and Abramovicz, 1996). 

6. Conclusions 

We conclude that plain ADAF solutions that don't respect the full set of sonic point 
conditions, including derivatives at sonic point, are incorrect. However, since ADAF models 
seem to offer explanation of relevant phenomena, to account for their existence, and in 
particular for the formation of a hot comptonizing corona, we propose that a subkeplerian 
thin flow is existing by itself, superposed to the standard keplerian disk. In accretion 
occurring in Active Galactic Nuclei this subkeplerian flows could easily come from the star 
cluster, that is known to have a very slow average rotation (Ho 1999). In low mass binary 
system this subkepler flow could come from a wind induced by the high energy radiation 
impinging on the normal star surface, or from gas recirculating in the binary potential well 
(Bisikalo et al. 1998). We are also testing the hypothesis that a physically motivated change 
in Z of the viscosity prescription can lead to differentially rotating flows. 
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Fig. 1. — Angular momentum per unit mass versus radial distance for a = .005, .01, .02. 
Fig. 2. — Mach number and sound speed versus radial distance for a = .005, .01, .02. 
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